la i y jaRS i Ti op illoois 

DIGIOM. COMPUTER 



LIBRAHT BOUTHE 2 •* 



TITLE 

TYPE 

lUKBER OF WORDS 

TEMPORARY STORAGE 

ACCURACY 

DURATIOH 



Poles 011*8 Sq^tloE, Li«1UDa9»«frtt|il^ NrKhod (ix)i Osily) 

Open 

Locations 6-6k, Grid points and key vorda may fill 

meicory « 

0-8 

10-5 

k PQ milliseconds, vhere P = number of Interior points 

Q <= mmber of tinea over grid 
For rectangular regions Q is roughly approxliBated by 



DESCRIPnOR 



5/login K* 



10 



where K* = ^o*^q • 1> 



OIq = [2//i^•(cos jc/p + cos x/q,) ]/[cos jt/p + cos it/q.] 

p + I =1 = number of columns, 

q^ + 1 » M = number of rows . 

This code gives the solution Y(x,y), corresponding to " 



given boundary values of Y and given interior values of p, of the 

2 
difference equation V- + Y- = /'h vhich approximates PbissoE's 
^ XX yy ' 

equation 



A 



7? 




' ^(x,y). 



The solution is found for all interior grid points of an arbitrary 
closed region with a square mesh of side h. The process used is that 
of Llebuiann and Frankel, which requires an initial approximation (la 
this case O) Y^ for each interior point (r,s) and then improves this 
approxination by the following iteration formula: 



-2- 
r^s r,s r-l,s r+l,s t^m^I r,s+l r,s 

where r = l,2,...,p - Ij s = 1,2,..., q - 1. For rectangular regions 

o{ - a = 2 -/yk ' [ cos </p + cos 3c/q] 
[cos jf/p + cos jt/q] 

and cX^ is computed by the "DetermineoC" routine, which also 
places it in location 5- This will suffice also for near-rectangular 
regions. For more peculiar regions omit "Determine oc " routine and 
set H(5) = lA "before " Evaluate Interior Points", 



Iteration stops when . , 

S I V"*^ - V° I < 5 X 10-^. 
r,s ^'^ ^'^ - 
The maximun] number of boundary and interior points depends on the 

complexity of the boundary but is limited to about 950 for very 

simple -boundaries. First the boundary should be enclosed in a 

rectangular grid of M rows and N columns, preferably with N > M. 

Known boundary points, of course, should correspond to points of 

the grid. If a given interior point has a sequence position i in 

the M X N grid ( reading from left to right and top to bottom ) then 

its neighbors have the position i - N, i - 1, i + 1,| i + N. Since 

the exterior points will not be stored it is necessary to find the 

actual addresses of the points in sequential position i, i-N, and i+N. 

This is done as follows: in preparing the problem a count is made of 

the number of points, with the sequences arranged in the same oilier as 

they are on the grid. These are called temporary key words and 

are placed on the tape after 00 65K as 20 nF OOF, 30 nF OOF, kO nF OOF 

for interior, boundary, exterior sequences, respectively, vfhere n is 

the number of points in a sequence. The total number of boundary 

and interior points plus the number of sequences of interior points 
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plus the number of sequences of interior points njay not exceed 9^1. 

The tape of temporary key words must be terminated with OOF OOF. It 

is placed in the memory in locations (65ff ) later occupied by the 

2 
function values and the ^h 's at the grid points t Temporary key 

words are operated on to form the addresses of grid points i» i*N« 

i+N, for each i which starts a sequence of interior points, and to 

place them into one word. These words are called the key words and 

are put into the memory in sequence starting at any location following 

the space required for the grid points. A suggested position is at 

address j = [99^ minus number of sequences of interior points]. Key 

-9 -.19 ..59 

words are of the form [mx2+px2 +qx2 ], where N(m), 

N(p), N(q) are the (i - N)-th, i-th, (i + N)-th numbers of the grid, 

respectively. 

All boundary values and P h 's are required to lie in the 

range - l/5 < V, Ph < I/5. The V*s anjaj)h 's are written and punched 

on the input tape as K(+) or S(-) followed by up to six decimal digits 

and are put in sequential order except that the character J must precede 

each boundary sequence (i.e. V's) and K must precede each interior 

sequence (i.e. /)h 's). The character F must be at the end of the input 

2 
data. Function 'values V and theP h 's are stored by the "Input" routine 

in the same word as follows: The sign digit distinguishes boundary points 
(1) and interior points (O); V occupies the next 20 binary digits; 
p h 's comprises the last 19 digits. Hence only one memory location 
is needed for each grid point, although the accuracy has been corres- 
pondingly reduced . 

2 
Printing of the solution is done in sequence, withioh 

for each interior point on the same line as the V for that point. 

"Input" parameters required are (l) N, the number of 

columns, to be placed in location 3 as N x 2* , and (2) j, the address 

•59 
for the start of the key words, to be placed in location. 4 as J x 2 . 

"Determine oC" , parameters are (l) p = N - 1 in location 3 a^ p x 2" , 



A- 



and (2) q = M - 1 in location if as q x 2 •. The "Print" parameter 
is s, the address of the last point to he printed, in location ^ as 



s X 2 
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The first grid point is at location 56. 
The conrplete routine requires the following tapes? 

1. Decimal Order Input (Xl) 

2. Temporary key words (in 65ff ) followed by OOF OOF 
and "Input" parameters. 

3« Routine to form key words (in 10-6^4-), followed by 26 
ION. 

k. Input routine (in 10 -^l-l), followed by 26 ION. 

2 

5. Boundary points V and Ph. 's, followed by F 

6. "Determine o( " parameters and routine (in 5-^6), 
followed by 22 l^^N. 

7. Input parameters again and routine for evaluating interior 
point (in 6-5^), followed \)y 26 6n. 

8. Print parameter, print routine (in 10-25), and Print 
routine (P2) (in 26 A5), followed by 26 ION. 

The order listed must be adhered to although only 2,5 and 
parameters need to be constructed for each problem. For non-rectangular 
regions, instead of (6) set N(5) = lA. Code (k) overwrites (3), (6) over- 
writes (h) and (2), (7) overwrites (6), and (8) overwrites (7). 

An example follows: 

N = 7, P = 5, 



M = 4, q = 3 



1 i 


3 2 


7 6k 12^^ 




12 


iS 


2k 


-5° 1 




18 


2k 


■ 50 


56 













"28 35 54 



91 



126 



152 



all values times 
h = .01 



For bouxidary points the given numbers are function values V; for interior 

2 
points ^>o h 's are given. All these numbers are in range. The temporary 

keys are 30 JF OOF, 20 kF OOF, 30 2F OOF, 20 kF OOF, 30 7F OOF, OOF OOF. 

' 2 
The list of boundary values and/? h 's is J, K, KDOOl, KOOO8, K0027, 

K006i+, KOI25, KDOOl, N, K0012, KOOI8, K0024, KDO3, J, K0126, KDOO8, N, KDOI8, 

K0024, KOO3, KOO36, J, KDI33, KOO27, KOO28, KOO35, KOO54, KOO9I, KD152,F. 
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Input parameters are N(5) «= 6 x 2 ■ and isik) = [998-2] x 2"^^ 

= 996 X 2" . "Determine oL" routine parameters are N(3) = 5 x 2"^^ 

and Hik) a 5 X 2'^^, The print parameter is k(^) = 79 x 2'"^^. 



RT: VlO/59 
DAT E September I3, 19^4 

PROGRAMMED B Y Richard F« King 



^ 



Igr 



LOCATIOH 


QiOIER 




ISCfTSS PAGE 1 




7(»M XET WORDS (POISS(S) 00 IQK 10»61|>F 





41 71* 








1^ 8r 




H(T) -W(8) -.0 


1 


50 (65 )F 

01 3F 


Ftem k6' 


4 " Curroct temporary key 


2 


LO 52L 
56 in 






3 


If 1 F 


By Vl 


Insert at cad of keys 




2lf 999P 




Read In "Inpwt" 


k 


LO 52L 


Frcam 2» 






LO 52L 




Test far interior seq.uence 


5 


56 IfTL 

t 


FO 






Ifl 5P 


Tes 


H(5) « N(6) =0 


6 


ko 6t 
L5 77 






7 


LO Of 




H(0) = [ninnber of interior points -f 




26 8L 




nwriber of boundary points 1 x 2 


8 


ifO F 
L5 IL 






9 


kS 26l 




Plant address of ctcrrent temporary 




lf6 kZL 




key 


10 


ID 51L 




Plant address of previous temporarj 




i^6 IIL 




key 


11 


50 F 


By 10M7 






01 5F 


From 17* 


Test for exterior sequence 


.12 


LO 55L 








56 18l 


Yes 




13 


11 3F 


No 




Ik 


L5 5F 
Sk F 
i+0 5F 




H(5) = [number of points we have 
backed v^] x 2'^'^ 


15 


LO 5F 








52 20L 


Yes 


Have ve backed up a row? 


16 


L5 til 


Ho 


Plant address of previous temporary 




LO 51L 




key. 



G2 



LOCMION , 


ORDER 




NOTES PAGE 2 G 2 


17 I 


i|6 IIL 
26 IIL 






18 


11 3F 
L5 6f 


from 12' 


N(6) = [ni:im"ber of exterior points "we 
have backed up] x 2" -^ 


19 


Sij- F 

4o 6f 






20 


22 13L 








L5 5^L 


From 15' 


A = [address of first mesh word] x 2" '^^ 


21 


L^ F 






22 ' 


i|0 IF 
LO 3F 
L^ 6f 




N(l) m X 2"-^'^ + p X 2"-^^ 


23 


00 lOF 








Lk IF 




y 


2k 


kO IF 

kl 5F 






25 


4i 6f 




n(5) = N(6) = 


26 


26 26L 
50 F 
01 3F 


By 9 
From 32' 


Test for exterior sequence 


27 


LO 53L 








36 33L 


Yes 




28 


11 "3F 


No 




29 


L5 5F 
S^ F 
kO 5F 


From 35 


N(5) = [number of points we have pro- 

1 -19 
gressedj x 2 


30 


LO 3F 




Have we progressed a row? 




32 35L 


Yes 




31 


L5 26L 


No 






lA 5IL 




Plant address of next temporary 


32 


k6 26L 
26 26L 




keye 


33 


11 3F 
L5 6f 


From 27' 


N(6) = I number of exterior points 
we have progressed] x 2 


3^ 


Sij- F 




, 



LOCATION 


ORDER 
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4o 6f 






35 


22 28L 






36 


L5 IF 
l4 3F 
LO 6f 


From 30' 


N(l) = m X 2"^ + p X 2"-^^ + q X 2"^^ 


37 


10 20F 
42 IF 




= key word 


38 


. L5 IF 








26 39L 




Store key word 


39 


40 S4 
L5 39L 


By 40' 




40 


l4 51L 
46 59L 




Plant addresses for next key word 


ifl 


46 3L 
26 42L 






42 


50 F 
00 3F 


By 9' 


Q = [number of points in present 
interior sequence ] x 2 


43 


11 3F 


From kQ 






L5 7F 


From 50' 


N(7) = [number of interior points 


44 


S4 F 




+ number of exterior points + 




40 7F 




number of boundary points] x 2 


^5 


L5 IL 








l4 51L 




Set address for next temporary key 


46 


46 IL 
26 IL 






i^7 


LO 52L 


From 5 


Test for boundary sequence 




32 48L 


NO 




48 


26 43L 


Yes 






11 3F 


From 47' 




49 


L5 8f 
S4 F 




N(8) = [number of e"j(ijerior points] 


50 


40 8f 
22 43L 






51 


00 IF 
00 F 




2-19 



LOCATION 


ORDER 
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52 


00 F 








00 IF 




2-59 


55 


00 F 








00 4f 




!♦ X 2-59 


^h 


00 56F 








00 F 




56 X 2-15 




26 ION 






INPUT 


, 00 lOK 


10~4lF 


F ■= End of input 
N - Interior sequence 
J - Boundary Sequence 
K - Positive Number 
S - Negative Number 





22 I5L 


From 18' 


Read in next character 


1 


75 IF 
00 39F 


From 7 




2 


^0 IF 
50 28L 






3 


L5 2F 
7^ F 






k 


S5 F 








40 F 


From 21,22' Code N3 j 


5 


81 4f 

50 28L 






6 


40 2F 
LO 26L 






7 


56 IL 
L5 F 






8 


66 IF 








L5 2k-L 




Test for interior point 


9 


56 25L 


Yes 






L5 29L 


NO 




10 


^h F 








10 IF 




Round off and position boundary 


11 


- 3S' 12L 




point number 



LOCATION 



ORDER 



NOTES 



PAGE ^ G 2 



12 

15 
Ik 

15 
l6 

17 
18 

19 
20 

21 
22 

23 
2i^ 

25 
26 

27 
28 



26 13L 
00 F 
ih 26l 
4o (56 )F 
L5 13L 

l4 51L 

k6 3IL 
L5 2F 
LO 27L 
3k 999F 
00 58F 
56 1% 
LI 2i^L 
i^-0 2ii-L 
26 L 
00 IF 
32 21L 
LI 3OL 
kO IF 
23 ^L 
L5 5OL 
kO IF 
23 ^L 
S5 F 
IM- 3IL 
10 21F 
36 I3L 
Ih 28L 
26 13L 
80 F 
00 lOF 
00 F 
00 i4f 

00 F 
00 lOF 



From 2ij-'25 
By Ik' 11' 



From 0' 
Yes 

Yes 



From 17 



From 19' 



From 9 



I Store point values and step 
address of store order. 



Test for F(end) 

Test for K or S 

Set N(23) negative or positive 

Start next sequence 
Test for K 

N(1) = »2"°^^ for ^ 



Form number 

N(i) = 2"^^ for K 



Form number 



Roundoff and position 

if 
(interior point number) x li' 



"1 + 10 X 2 
Ik X 2°^^ 
10 X 2"^^ 



"39 



LocAnoaf 


ORDER 


NOTKB 


PAGE 6 


29 


00 F 
8o F 




2-20 




30 


00 F 
00 IF 




2-59 




51 

1 


00 IF 
00 F 
26 ION 




2-^9 




N(5) ' 


P X 2"^^ 








K(i^) « 


4 X 2-59 




Detemdne ^ 




00 5K 


5-15F 











51 lOL 
00 2F 




-39 - 
A Q = « X 2 ^^ 




1 


22 IL 
66 3F 




A - x/p 




2 


85 F 
50 2L 


1 


H(3) - (coa it/p) 1/2 




5 


26 i6f 
L5 2r 








k 


if 3r 

51 lOL 








5 


00 2r 










66 kT 




A - JcA 




6 


85 F 

50 6l 








T 


26 i6f 
L5 2F 




A - (coa n/v) I/2 




6 


lA 35F 
ifO 5F 




N(3) - (coa jt/p + coa 


«/4) 1/2 


9 


26 999F 




End of interlude 






50 3F 






^ 10 


kOF 00 2855 








9816 3397 J 


x/l^ 






00 i6k 






SiB«, 


CoilM (Tl) 16 - if6F 
26.J3L- — _ 







LOCATION 



ORDER 



NOTES 



PAGE 



G^2 



0015K 15-^2p 
7J 5F 

1^0 5F 
L5 7L 
LO 3F 
10 2F 
50 2L 
26 25F 
49 IF 
il- LO 2F 

66 3F 

5 S5 F 
40 5F 

6 2^ 999E 
OF F 

7 7L i^095F 
LL 1|095F 

00 23K 23 " 52F 

22 l^N 
Evaluate Interior Points (Poisson) 



00 6k 




6" 54F 

41 5F 
L3 S4 
32 35L 
L5 Si| 
40 2F 
h6 7L 
1^6 30L 
h2 i6l 
L4 41L 
46 15L 
LO 42L 
46 i8l 
10 30F 

42 19L 



N(5) _ /cos ^/p + cos rt/q j 



= l/l^. [1 „ rcos it/p + c os jt/qV J 






(5) =^ 




By 9\59 

By 10, 39' 
From 35 



[ (cos gt/p + cos :t/q) l/2 ]' 
Read in next routine 
Stop 

1 - 2-39 

Square Root Routine (R1) 

N(5) =<^ 

N(3) = 

Test for end of grid 

N(2) = current key word 



Plant p,q 
Plant p+1 
Plant p-1 

Plant m 



LOCATION 
7 

8 

9 

10 
11 
12 

14 

15 
16 

17 
18 

19 
20 
21 
22 

23 
2i^ 



ORDER 
L5 (p)F 
i^O hY 
36 IIL 
L5 L 

l4 4il 

42 L 
42 IL 
22 L 
50 46l 
00 21F 
10 3F 
50 4f 
JO 47L 

54 F 
00 2F 
40 F 

50 (p+l)F 

JO 47L 

55 F 

50 (q)F 

JO 47L 

S4 F 

50 (p»l)F 

JO 47L 

S4 F 
50 (m)F 
JO 47L 

S4 F 

LO F 

00 IF 

40 F 

50 F 

75 5F 
L4 45L 
4o F 



NOTES 



PAGE 8 2 



By 2' 



Yes 



From 



By 4 



By y 



By 5' 



By 6 



N(iSf) A = current function value 
Test for boundary point 

Plant key -word addresses 



End o£ interior sequence 



A = 1/2 (/h^A) 



A = N(0) = 1/2 (+4 0^ + f^) 



A = 1/2 (4 0"^) 



n -i.^n ,7,2 



N(P) = 40"^ = :40''» / h' 



N(0) =c?r[40"''-40^ y^h^] 



LOCATION 


ORDER 
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10 IF 




r -} 


25 


50 ij-F 
JO i|-7L 




N(l) = 1/2 |0° +0{[ll]j = 1/2 0°+^ 


26 


Sk F 
i^-0 IF 






27 


50 4f 

JO ij-6L 






28 


S5 F 
50 IF 




H(p) =1/2 [0°+^ + 2-2°x/Ph2] 


29 


JO 47L 

S4 F 






50 


kO (p)F 

50 F 


By 3 




31 


JO hQL 
S7 F 




N(3) =E H^+j_(p) - H.(p) 


32 


lA 3F 
kO 3F 






33 


36 34l 




Test for overflow [n(43) < O] 




kO i^-3L 


Yes 




5^ 


L5 2F 


From 33 






I^ i^■lL 




Step m^ p, q by 1 


35 


26 2L 








LI i^3L 


From 1 


Test for overflow 


36 


36 38L 








L5 ii-l+L 


NO 




57 


LO 3F 
3k 999F 




Test for end 


38 


i^9 i]-3L 
L5 i^-OL 


From 36 


Reset N(il-3) > 


39 


k2 L 




Reset addresses and start over on 




il-2 IL 




grid 


ifO 


26 L 








OF Sk 




Address used 


41 


00 1025F 

00 1025F 




2-9 + 2-19 ^ 2-29 ^ 2-39 



LOCATION' 


ORDER 
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k2 


00 1026F 
00 1025F 




2-9 + 2 X 2-15 + 2-29 + 2-59 


^^ 


kO F 








00 F 




1/2 


hk 


OOF 00 0005 
0000 0000 J 


5 X 10 ^ 


^5 


00 F 
80 F 




2^20 


hG 


00 F 








7L il-095F 




00 »- 


001 — 1 


hi 


7L i^095F 










80 F 




ail"- 1 


100( ^^ 


kQ 


LL i^095F 










00 F 




11 -.« 1 00-— 




26 6n 






Print (] 


Poisson)00 lOK 




10=>25F 





92 129F 


From 11 






50 (56 )F 


By 9' 




1 


JO li^-L 




Extract, position + print function 




00 il-OF 




{ 


2 


50 5F 
50 2L 




values 


5 


26 26F 








L5 (56 )F- 


By 10 




k 


36 5L 


m 


Test for boundary point 




22 8L 


Yes 




5 


22 5L 


From ^ 






50 (56 )F 


By 10' 




6 


JO 15L 
00 6OF ' 




■ 

Extract, position, and print /^h 's 


7 


50 3F 
50 7L 






8 


26 26F 








19 38F 


From k ' 




9 


lA L 




Set addresses for next point 



LOCATION 


ORDER 
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k2 L 






10 


k2 3L 
i]-2 5L 






11 


LO 15L 




Test for end 




52 12L 


Yes 




12 


26 L 


NO 


Print values for next point 




OF F 




Final stop 


13 


92 129F 








50 ISif 




End constant 


Ik 


7L 4095F 








80 F 




Oil -- 1 


100 — 


15 


00 F 










7L ^095F 




000 — - 


-11 — 1 




00 26k 








Print Routine (P2) 






26 ION 







